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^ • Abstract 

o 

f-j I We expand our study on the dispersion of polarization angles in molecular clouds. 

^ ■ We show how the effect of signal integration through the thickness of the cloud as 

well as across the area subtended by the telescope beam inherent to dust continuum 



(N 



measurements can be incorporated in our analysis to correctly account for its effect on the 
measured angular dispersion and inferred turbulent to large-scale magnetic field strength 



> 

04 ■ ratio. We further show how to evaluate the turbulent magnetic field correlation scale 

from polarization data of sufficient spatial resolution and high enough spatial sampling 
rate. We apply our results to the molecular cloud OMC-1, where we find a turbulent 



O ' correlation length of (5 « 16 mpc, a turbulent to large-scale magnetic field strength 

^ ■ ratio of approximately 0.5, and a plane-of-the-sky large-scale magnetic field strength of 

approximately 760 /iG. 

■ 

•i-H - 
V . 
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1. Introduction 

The observational determination of the turbulent energy content within the magnetic field is 
important for understanding the role of magnetic fields in the star formation process. It provides, 
for example, some measure of the amount of turbulent energy contained in the gas while the further 
determination of the turbulent to large-scale magnetic field strength ratio can in principle be used 
to evaluate the large-scale ma gnetic strength with the so-called Chandrasekhar-Fermi equation 



|Chandrasekhar fc Fermilll953l ). 
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The determination of the turbulent to large -scale magnetic field strength ratio was the subject 
of a recent publication (IHildebrand et al.l l2009l : hereafter Paper I) where it was shown how this 
parameter can be precisely extracted from dust continuum polarization data through a careful 
analysis of polarization angle differences as a function of the distance between pair of points where 
measurements were made (i.e., the angular dispersion function). It was described how the evaluation 
of the turbulent polarization angular dispersion can be achieved without assuming any model for 
the large-scale magnetic field component about which this dispersion is calculated. This is an 
important development since no model will perfectly fit the true morphology of the large-scale 
magnetic field. Fits to a model will therefore lead to inaccurate estimates of the angular dispersion. 
This error would then be propagated in the determination, for example, of the large-scale magnetic 
field strength when the Chandrasekhar-Fermi equation is used. 

In this second paper on the subject we generalize our analysis of Paper I by including the 
process of signal integration through the thickness of the cloud as well as across the area subtended 
by the telescope beam inherent to dust continuum measurements. It has long been recognized that 
the amo unt of angular dispersion measured in a polarization map is reduced by any integration 



process (IMvers fc Goodman 



1991 



cal simulations (Qstriker et al 



as has since been studied and demonstrated through numeri- 
200l|; IPadoan et aklbonil : iHeitsch et aklbonil : Ixudoh fc Basulbond : 



Falceta-Goncalves et al.ll2008ll. The effe ct of integration through the thickness of a cloud has been 



considered by lMyers fc GoodmanI (Il99ll ) in their studies of the optical polarization of dark clouds, 
while the further inclusion of integration across the area subtended by the tele scope beam and ensu - 
ing consequences on me asure ments were investigated throug h simulations by lHeitsch et al.l (|20Qll ). 



Wiebe fc WatsonI (120041 ). and iFalceta-Goncalves et al.l (120081 ). 



We will start with a generalization of the problem considered in Paper I by deriving the cloud- 
and beam-integrated dispersion function in Section [2l which will then be solved for the special 



case of Gaussian turbulent autocorrelation and beam profile functions. In Section [3l we a 



analysis to dust continuum polarization data obtained with SHARP |Novak et al 



jply our 



2004; 



Li et al 



20061 . 12OO8I ) for the molecular cloud OMC-1. We show how the turbulent correlation length scale for 



this cloud can be evaluated with the corresponding polarization data. We then use these results to 
determine the number of independent turbulent cells contained in the column of dust probed with 
our measurements and calculate the turbulent to large-scale magnetic field strength ratio corrected 
for the signal integration process. We provide a detailed discussion of our results in Section H] and 
end with a summary in Section [H Detailed derivations resulting in the relations and functions used 
in these sections, description of the data analysis, as well as a list of variables and functions will be 
found in the appendices at the end of the paper. 
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2. Analysis 

2.1. The Cloud- and Beam-integrated Angular Dispersion Function 

In Paper I the following equation was introduced (see their Equations [A4] and [A20]) for the 
analytical derivation of the dispersion in polarization angles within a turbulent molecular cloud 



[(BMx»(-BMx + 0)l'"' 

where A$ {£.) = $ (x) — $ (x + £) is the difference in the polarization angle $ measured at two 
positions separated by a distance £ and (• • • ) denotes an average. As was then shown, Equation 
((T|) (and others that derive from it) applies equally well to a three-dimensional magnetic field or its 
two-dimensional projection onto a plane; for the purpose of this paper we are considering B (x) to be 
the plane-of-the-sky projected magnetic field, as usually probed with dust continuum polarization 
measurements. The distance £ is also confined to the plane-of-the-sky, unless otherwise noted. We 
assume that the magnetic field B (x) is composed of a large-scale, structured field, Bo(x), and a 
turbulent (or random) component, Bt (x), such that 



B(x) = Bo(x)+Bt (x). 



(2) 



We must note, however, that the magnetic field direction inferred from polarization data is 
actually the result of some averaging process as one integrates through the thickness of the cloud 
along the line-of-sight as well as across the area subtended by the telescope beam. We therefore 
first define the normalized magnetic field vector 



b(r,z) 



1/2 



(3) 



and infer a mean direction for the cloud-integrated magnetic field through the following weighted 
integral of b (r, z) 



h(r) 



jjH{r-s.) ^J%{R,z)h{sL,z)dz 



(4) 



The two-dimensional (convolution) integral in Equation ^ is over all space such that for the three- 
dimensional position vector x, r is the two-dimensional polar radius vector on the plane-of-the-sky 
and z the depth within the cloud. That is, 



X = re^ -|- zGz 



(5) 
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with and the unit basis vectors along r and the z-axis (which is oriented along the line-of- 
sight), respectively. The beam profile density is denoted by H{r), while the weighting function 
F (r, -z) > is the polarized emission associated with the magnetic field b (r, zj^. Please note that 
even though b (r, z) is normalized b (r) is not, and while the normalization by A in Equation ^ 
is not essential it is included for convenience. Moreover, the quantity A is for the maximum depth 
of the cloud along any line-of-sight; the detailed behavior of the F (r, z) function thus ensures that 
Equation ^ is exact even when the depth of the cloud is expected to vary with position on the 
plane-of-the-sky. 

The normalization of the magnetic field vector through the cloud is warranted because the 
amount of polarized emission in a given region is not a function of the strength of the magnetic 
field itself. Because of this we must now consider a slightly different relation for determining the 
dispersion function. That is, we replace Equation ([T]) with 

(cos [A^ m ^ ; \' ' (6) 

(f (r)) (f {r + £))] 

In what follows we introduce a level of idealization necessary for obtaining a quantitative mea- 
sure of the turbulent component of the magnetic field in molecular clouds. We assume stationarity, 
homogeneity and isotropy in the magnetic field strength, as well as statistical independence between 
its large-scale and turbulent components. We therefore have the following averages at points x and 

y 



(Bo(x)) = Bo(x) 

(Bt(x)) = 

(Bo (x) • Bt (y)) = (Bo (x)) • (Bt (y)) = 0, (7) 

and 



(B^(x)> = {BUy)) = {Bl) 

(B?(x)> = {BUy)) = {B^). (8) 

It is straightforward to show that the homogeneity in the field strength renders the field nor- 
malization of Equation ^ inconsequential in Equation ([6]) and we can therefore equally write 



^The measured linear polarization orientation is normal to that of the associated plane-of-the-sky magnetic field 
when detecting dust continuum emission (at submillimeter wavelengths, for example), not parallel to it as could be 
inferred from Equation ([4|). This is irrelevant to our analysis, however, as we are using the polarized emission as a 
weighting function in defining a mean orientation for the integrated magnetic field. 
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(cos [A$ {C)]) 



(B(r)-B(r + . 



(r)\/s' {v + l) 



1/2-- 



(9) 



with the cloud- and beam-integrated magnetic field 



B(r) 



H(r-a) 



A 



F (a, z) B (a, z) dz 



(10) 



and where the assumed isotropy in the distance {I = |£|) was incorporated. Equation ([9]) is the 
relation we will use to estimate the angular dispersion function. 



2.2. The Integrated Magnetic Field Autocorrelation Function 

In view of the assumed stationarity and isotropy the integrated magnetic field autocorrelation 
function <^B-B(^)) = (B(r) •B(r + ()) can be expressed as (see Eqs. |Al| - [A5j in Appendix [A|l 



(B-B(^)) 




H (a) H (a' + i) 



^ f (l - ^) {v, u) du 

J 



d a d a, 



fill 



where we introduced the (assumed statistically independent) autocorrelation functions for the (non- 
integrated) magnetic field and the polarized emission 



Rb{v,u) = (B(a,z) •B(a',z')> (12) 
Rf{v,u) = {F {a, z) F {a' , z')) (13) 

with u = \z' — z\ and f = |a' — a|. These autocorrelations can be further broken down through the 
decomposition of the magnetic field and the polarized emission into their respective large-scale and 
turbulent components (see Equation [2]) 



B(a,z) = Bo (a, z) + Bt (a, z) (14) 
F(a,z) = Fo (a, z) + Ft (a, z) . (15) 

Upon assuming the same statistical properties for the components of F (a, z) as those expressed in 
Equations ([7]) for B (r, z), Equations (fT2]l and (fT3]) transform to 
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Rf (v, u) = Rf,o {v, u) + RF,t (v, u) 



(16) 
(17) 



with 



RB,jiv,u) = (Bj(a,z)-Bj(a',z')> (18) 
Rfj {v, u) = {Fj (a, z) Fj (a', z')) (19) 

where j = '0' or 't' for the large-scale and turbulent components, respectively. 

The solution to our problem is reduced to solving Equation ([TT]) given the different autocorre- 
lation functions; the equation to be used for determining the dispersion function then becomes 

(cos [A^ {£)]) = , (20) 

^ ^ ^ (B-B(0)> ^ ' 



2.3. Solution Using Gaussian Turbulent Autocorrelation and Beam Profile Functions 

In order to solve Equation ((20]) we must specify, at least to some level, the characteristics of 
the different autocorrelation functions, as well as the telescope beam profile. We therefore define 
the following 



Rb{v,u) = RB,o{v,u) + {B^)e-(-"+'-")/''" (21) 
Rf{v,u) = RF,o{v,u) + {F^)e-(-'+-')/^''\ (22) 

where the magnetic field and polarized emission autocorrelations each have large-scale and turbulent 
(second term on the right-hand sides) components. The correlation length scales for the turbulent 
magnetic field and polarized emission are 5 and 6' , respectively, and, as is implied through Equations 
([2T]) and (f22l) . we assume the turbulence to be isotropic. Furthermore, 6 and 6' are taken to be 
much smaller than the thickness of the cloud (i.e., 5 <C A and 5' <C A). The beam profile is defined 
with 



where W is the beam "radius". 
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Under these constraints and definitions Equation (fTT]l can be analytically solved to yield (see 
Eqs. [A6]-[A8] Appendix El) 



(B.B(£)> ^ {Bi){FS)i{aii)) + V2^- 



6^ 



(52 + A_ 



-^2/2(52+21^2 



y/3 



((5"2 + 2W^) A 



((^'2 + 2T^2) 



g-£2/2(5"2+2W2)| 
"' g--^2/2(y2_^2VK2) 



(24) 



with 



and the normalized large-scale function 



Rb,o (0,0) 
Rf,o (0,0) 
66' 



(25) 
(26) 

(27) 



{a {£)) 



i?(a)i/(a' + £)||^' 




It is expected that this large-scale function will usually dominate the other terms in Equation 
([241) . since these all result from the averaging of the turbulent components through the column of 
dust probed by the corresponding polarization measurements. Indeed, it is apparent that each of 
the turbulent contributions in Equation scales with a term such as 



N 



-1 



27r5^ 



(52 + 2Ty2) A' 



(29) 



where N is nothing more than the number of independent turbulent cells (for the magnetic field in 
this case) contained in the column of dust probed observationally, as could have been intuitively 
guessed. We also note that for cases where the telescope beam radius were to be much smaller than 
the size of a turbulent cell (e.g., 6 ^> W) we recover 



N 



27r6 



(30) 
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The number of independent turbulent cells would thus be accounted for by those cells that lie along 
the line-of-sight thro ugh the thickness of the c loud at a given point on its surface, as would also be 



expected intuitively (IMyers fc Goodmanlll991 



Inserting Equation (p4|l through (|28]) into Equation (1201) we can write 



1 - (cos {£)]) 



(B-B(O)) 



[{am) -{aim 



5^ 



+ 



((5"2 + 2Ty2) A 



1 _ g-^V2(52+2W2y 
1 _ g-^V2(<5"^+2l^2y 



+\/2^' 



(5'2 + 21^2) A 



1 - e 



-^V2('5'^+2VK2 



(31) 



The first term within parentheses on the right-hand side is, as was previously stated, due to the 
large-scale structure in the magnetic field and the polarized emission; it does not involve turbulence. 
We can expand this term using a Taylor series with 



oo 

(a(0))-(Q(£)) = ^a2/2^ (32) 

where the summation is performed only on even values for 2j, since {a{tj) is isotropic in ^ (see 
Section [Q in Appendix (A]). 

Although the remaining terms in Equation ([3T]) are all due to turbulence, their respective 
contributions have, in part, different origins. More importantly, it should be clear that the total 
turbulent component is not only due to turbulence in the magnetic field but can also arise from the 
presence of turbulence in the polarized emission. For example, the third term within parentheses 
on the right-hand side of Equation (fSTI) can be interpreted as a contribution due to random changes 
in polarized emission at different positions in the cloud where the magnetic field also changes 
orientation. It follows that the measured (i.e., integrated) orientation associated with the polarized 
emission, and therefore the deduced magnetic field orientation (i.e., the polarization angle), will 
also accordingly fluctuate randomly. 

The last term in Equation (j3ip seems to imply that the turbulent component in polarized 
emission will contribute to the turbulent angular dispersion even in cases where the magnetic field is 
uniform (and does not have a turbulent component). But this is an artifact of the way we analytically 
evaluate the dispersion function through Equation (j3ip . which does not perfectly mimic the way 
polarization measurements are accomplished (see Section H]). In reality, polarization measurements 
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are made on a point by point basis, and the dispersion function is calculated through an average 
of (the cosine of) angle differences as a function of the displacement i (see Eq. [1] in Paper I) not 
through the evaluation of an autocorrelation function such as given in Equation (fTTj) . It follows 
that a perfectly uniform magnetic field could never lead to a measurable angular dispersion; we will 
therefore not include in our analysis the corresponding contribution in Equation pTI) . 

Now consider the following 



(B-B(O)) 



{a{Q)) 



iJ(a)ff(a') 



Rf.o («,m) 



-Rb.q {v,u) 



du \ (fa'cPa, (33) 



where we have taken advantage of the fact that we expect that the large-scale component (a (0)) 
dominates the turbulent terms in Equation (f24l) (when ^ = 0), as was previously stated. Because the 
autocorrelations present in the integrand of the one-dimensional integral are normalized it follows 
that this integral, which we denote as A (|a' — a|) /2, will always be less than A/2 and thus 



(B-B(O)) 



Ip2\ 
I 




1 

A 



A' 



i7(a)F(a')^(|a'-a|)(i2a'd2Q 



(34) 



where the equality to unity only holds when the large-scale magnetic field and polarized emission 
are both uniform. The quantity A' can be interpreted as the effective thickness of the cloud; this 
will be discussed in more details in Section [321 

Taking these considerations into account, Equation ([3T]) becomes 



1 - (cos [A$ (€)\, ~ V27r 



53 



+ 



I 



(52 + 2H^2) ^ 

5"3 



(5"2 + 2W^2) ^ 



1 _ g-^V2(52+2VK 
1 _ g-^V2(<5"^+2H/2) 



+ E«2/'' (35) 



with a'gj- = (A/A')a2j. If we limit ourselves to small enough displacements such that I is less than 
a few times the beam radius W , then the large-scale term becomes small enough to be adequately 
described by the first term in the Taylor expansion of Equation ([32]) and Equation ([35]) can be 
approximated to 
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1 - (cos[A$(^)]) ~ V27r 



or alternatively 



[Bl) 



6^ 



(52 + 2W^) A' 
y/3 



1 _ g-^V2(5"^+2Vt'2) 



(36) 



(A$2 (£)) ~ 2\/2^' 



where m? = 2a 



2- 



'B' 



+ 



{Bl) 



6^ 



(52 + A/ 



(J//2 + 2VF2) 



1 _ g-^V2(52+2VK2 
1 _ g-<?V2(5"2+2iy2 



+ m 



2/)2 



(37) 



3. Results 
3.1. The Polarized Emission 

Ideally one would intend to use Equation ((36]) for the dispersion function (or alternatively Equa- 
tion [37]) to determine the ratio of the (square of the) turbulent to large-scale magnetic field strength 

) / (^o)- shown in Paper I, this dispersion function is readily evaluated from polarization 

data and the aforementioned ratio could be used to calculate, for e xample, the large-scale magneti c 
field strength (-Bq)^^^ through the Chandrasekhar-Fermi equation ( Chandrasekhar fc: Fermi 1953 ). 
It should be clear, however, that in order to precisely achieve such a goal we must find a way to 
determine other parameters such as the turbulent correlation scales 5 and 5', the effective cloud 
thickness A', and the ratio of the turbulent to large-scale polarized emission (-Ft^) / {Fq)- 

The polarized emission, in particular, is problematic. This is because we do not have direct 
information on the weighting function F (r, z) or its integrated counterpart 



F(r) 



1 /"^ 
a/o 



(38) 



Although the integrated polarized flux P (r) i s contained in pola rization da ta sets obtained 



with dust continuum polarimeters such as H ertz (IDotson et al.ll2009l ). SCUBA (jMatthews et al 



20091 ). and SHARP (jVaillancourt et al.ll2008l ). it does not correspond to the quantity defined in 
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Equation (f38]l . Instead the polarized flux is observationally determined through the measurement 
of the integrated Stokes parameters 



Q{r) 
U{r) 



H{r-a) 
H{r-a) 



A 



A 



Q (a, z) dz 



U (a, z) dz 



d^a 



d^a 



(39) 
(40) 



with 



P(r) 



{r) + U^ (r). 



(41) 



Although these parameters will also exhibit large-scale and turbulent components, they do not 
provide us with any means for disentangling the turbulent contributions (-Bt ) (i^t^) due to the 
magnetic field and the polarized emission they contain. We must then resort to some approxima- 
tion(s) if we are to make any progress. 

Perhaps the most obvious difference between the two turbulent terms on the right-hand side 
of Equation ([351) is that the first one is of first order in the square of the turbulent to large-scale 
magnetic field strength ratio while the other is of second order. If we assume that such a ratio is 
a fraction of unity, then the second order term can be neglected. This amounts to neglecting the 
contribution of the turbulent polarized emission in Equation (f35l) : this is the line of attack we will 
use from now on. We therefore write that 



1 - (cos [A$ {I)]) ~ V27r 



(52 + 2W^) A' 



00 

^_^-iV2{5^+2W-)^^J2a',/^, (42) 



while for displacements i less than a few times W we keep only the first i term in the Taylor 
expansion 



1 - (cos [A$ (£)]) ~ V2^T^ 



(<52 + 2W^) A' 



1 - e 



(43) 



or 



A$^ (^)) ~ 2V27r 



5^ 



(52 + 2W^) A' 



1 - e 



(44) 



Whether it is appropriate or not to neglect the turbulent polarized emission is an open question, 
which could perhaps be adequately investigated through simulations. 



3.2. The Effective Cloud Thickness 



We previously defined the effective cloud thickness A' through the following relation 



{am 




R-B,o iv,u) 
{Bl) 



A 



< 1. 



(45) 



We can get a sense of the nature of A' by expressing (a (0)) with the Fourier transform of (a [€)) 
at £ = (see Section [A. 21 in Appendix [A] for more details) 



^0 (k. 



■V 1 ) 



smc 



kuA 



(46) 



where TZq {ky , ku) is the Fourier transform of TZo{v,u) = RFfl{v,u) RBfiiv,u) and sine (x) = 
sin(x)/x. It should be apparent that Equation ([46]l includes several effects that set the value 
of A'. For the idealized case where the autocorrelation Tlo{v,u) is uniform across the cloud, its 
Fourier transform TZq (k^,, k^) is proportional to a Dirac distribution (more precisely, TZq (k„, /c„) = 
(27r)^ (-^0 ) (^o) ^ (kt;) ku)) and A' = A, as was previously stated. For more realistic cases, however, 
TZq {ky, ku) will have a finite width along ku and k„. As is also discussed in Section [AT2I the (square 
of) beam profile \\H (kt,)||^ will filter IZq (k^,, ku) along k^ and therefore reduce the value of (a (0)) 
(i.e., A' < A) through the exclusion of spectral modes located outside the bandwidth subtended by 
(k^)ll^. A more important filtering effect due to the finite spectral width of TZq (kt,,/c„) along 
ku is expected, however, as it is severely truncated by the integration process through the cloud 
thickness A. This is clearly assessed by the presence sinc^ (/ctjA/2) in Equation (l46ll and the subse- 
quent integration on ku- Because of the assumed large difference between the cloud thickness and 
the beam radius (i.e., 14^ <C A) we expect that this effect will dominate in determining the value of 
(«(0)>. 

If we now introduce Afe^ the spectral width of TZq (k^,, ku) along ku with the following definition 



2Afc„7^o (k„,0) = y"7^o(k 

V 1 ) dku 5 



(47) 



then we can further approximate (since Aku 3> A ^) 
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TZq (kt,, ku) sinc^ ' 



^ dku ~ 7^o(kt,,0) J sinc^ 



~ ^7^o(k^„0) 
~ a" 2AA;„ 



(48) 



Inserting this last relation into Equation (l46l) . with the assumption that the filtering due to the 
telescope beam is negligible compared to that due to integration through the thickness of the cloud, 
yields 



(a(0)) ~ — 




d-'k 



(49) 



We then find 



A/c„ ~ — . (50) 

Because of the way the Fourier transform relates the width of TZq (k^, ku) to that of TZq {v,u), A' 
can therefore be advantageously interpreted as the width of the large-scale autocorrelation function 
TZq {v,u) (if we assume isotropy, then the width is the same along u oi v). 

Unfortunately, we do not have access to TZo {v, u) when mapping the polarization of dust emis- 
sion in a molecular cloud. However, we can probably get a decent approximation for A' through the 
shape of the autocorrelation function of the aforementioned cloud- and beam-integrated polarized 
flux i-'(r). This function is defined with 



(p\l))^{P{v)P{v + l)). (51) 

Albeit this is an integrated quantity, it does provide a good sense of the proportion of the cloud that 
contains the bulk of the polarized flux, which is used as the weighting function for the magnetic field 
in our analysis (see Eq. [lO]). We also note that this autocorrelation is a function of the distance I on 
the surface of the cloud, not through its depth. But this is consistent with the isotropy assumption 
from which it is expected that a molecular cloud will have similar characteristics through its depth 
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with the 



and across its surface. It therefore seems reasonable to associate the width of 
effective depth of the cloud. 

For the case of OMC-1 we use previously published Hertz datsl^ ( Houde et al.ll2004 : Hildebrand et al 



2003 ) to evaluate the effective depth of the cloud. The result is presented in Figure [TJ where the 
normalized autocorrelation function of the integrated normalized flux is shown. As indicated on the 
graph, we have arbitrarily chosen the width at half magnitude as the value for A'; we will therefore 
set A' « 3'.5 for OMC-1 for the calculations that will follow. Also notably, we redefine N as follows 
(see Eq. [29]) 



N 



+ 2W^) A' 



(52) 



3.3. The Turbulent Correlation Length Scale, the Turbulent to Large-scale Magnetic 
Field Strength Ratio, and the Large-scale Magnetic Field Strength 

Having estimated the effective depth of the cloud A' from the autocorrelation function of the 
integrated normalized flux, we are now in a position to determine two fundamental parameters that 
characterize magnetic fields and turbulence in star-forming regions: the turbulent correlation length 
6 and the (square of the) turbulent to large-scale magnetic field strength ratio (-B^) / (Sq)- To do 
so we refer to Equation (I43p . which we write again here for convenience 



1 - (cos [A$ {£)]) 



53 



(<52 + 2W^) A' 



(53) 



which is valid when the displacement I is less than a few times W . 



Our plan consists of using the previously published 350-/im SHARP polarization map of OMC-1 
(jVaillancourt et al.ll2008l ) to evaluate the left-hand side of Equation ((531) and fit our solution to the 
problem (i.e., the right-hand side) to the data. There are only three quantities to be simultaneously 
fitted for: 5, (-B^^ / (-Bg), and a'g, the first two being the parameters we are most interested in at 
this time. More details concerning our data analysis will be found in Appendix IbI 

We show in Figure [2] the result of our non-linear fit to the aforementioned data. In the figure, 
the top graph shows the fit of Equation ([53]) (solid curve) to the data (symbols) when plotted as 
a function of . The broken curve does not contain the correlated part of the function, i.e., the 
function 



^We use data obtained with Hertz instead of SHARP as the former cover a larger spatial extent and allow us to 
determine (P {£)) over a large enough distance. 
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is displayed to better visualize the integrated turbulent contribution (i.e., the first term) to the 
dispersion in relation to the large-scale component (i.e., a2^^) when £ is less than a few times W. 
The middle graph of Figure [2] displays the same information as the top graph but plotted as a 
function of £. In the bottom graph we show the correlated turbulent component of the dispersion 
function (symbols) obtained by subtracting the data to the broken curve in the middle graph, which 
in our model (solid curve) corresponds to 



Finally, the broken curve shows what would be the expected contribution of the (assumed 
Gaussian) telescope beam alone to the width of the turbulent component (i.e, when <5 = in 
the argument of the exponential). Although we performed our analysis for a Gaussian turbulent 
correlation function, it is rather unlikely that the autocorrelation function of the turbulent magnetic 
field component fits this model. We therefore did not use the first few points (i.e., for £ < Of 2) to fit 
Equation (f53]l to the data and concentrated on larger values of £ where it is more likely to obtain 
a reasonable fit. Indeed, there is evidence from the first few points, when £ < 0^2, in the middle 
and bottom graphs of Figure [2] that the Gaussian turbulent autocorrelation function assumption is 
incorrect, as expected. 

A comparison of the expected contribution of the telescope beam to the correlated turbulent 
component of the dispersion function (broken curve in the bottom graph of Figure ^ with the 
accompanying data reveals the imprint of the finite turbulent correlation length scale 6 in the 
relative excess detected in the O'.l < ^ < 0^4 range. This imprint is also evident from the non-zero 
integrated ratio of the (square of the) turbulent to large-scale magnetic field strength ratio (i.e., 
6^ (0) from Equation [55]) that is evident in the data through the intercept of the broken curves at 
£ = in the top and middle graph of Figure [2l As is seen in Equation (f55l) this intercept will go to 
zero in the limit where (5 — > when the number of turbulent cells subtended by the telescope beam 
A'" tends to infinity (see Equation [52]), i.e., in cases where the turbulent component is basically 
completely integrated out. 

The results from our fit of Equation ([531) to the dispersion data are summarized in Table [H 
Most notably, we measure the turbulent correlation length to be (5 f« 16 mpc (or 7'.'3 at 450 pc, 
the distance we adopt for OMC-1), which implies that there are on average N 21 independent 
turbulent cells contained within the column of gas probed by our telescope beam. Furthermore, 




(55) 
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since our fit for the square of the "integrated" turbulent to large-scale magnetic strength ratio yielded 
6^ (0) « 0.013, then it follows that 



Nb"^ (0) 
0.28, 



(56) 



where we used Equations (|52]) and (l55l). This quantity can be inserted in the Chandrasekhar-Fermi 
equation (IChandrasekhar &: Fermilll953l ) to evaluate the strength plane-of-the-sky component of the 
large-scale magnetic field 



« 760 nG, 

where we used the same values as in Paper I for the mass density p (i.e., a gas density of 10^ 
cm~^ and a mean molecular weight of 2.3) and the one-dimensional velocity dispersion a (v) (i.e., 
1.85 km as obtained from a representative H^^CO J = 3 — > 2 spectrum). 

We are now in a better position to appreciate the importance of adequately taking into account 
the signal integration process. Indeed, failing to do so while using the Chandrasekhar-Fermi equation 
would imply multiplying the value obtained in Equation (f57l) by \/iV 4.6. This would yield a 
magnetic field strength of approximately 3.5 mG, a value that is certainly prohibitively high for 
OMC-1 (see Section |43|) . 



/ r2\ 



-1/2 



(57) 



4. Discussion 

4.1. The Turbulent Power Spectrum 

The determination of the turbulent correlation 5 length is important for the characterization 
of turbulence in molecular clouds. It is therefore desirable to compare our result of 5 16 mpc for 
OMC-1 with other independent techniques or analyses that seek to evaluate this quantity, or others 
related to it, either observationally or theoretically. 

A parameter that is very closely related to 5 is the ambipolar diffusion scale 5ad at which the 
ionized and neutral components of the gas decouple, and that may determine the cutoff wavelength 
of the power spectrum of the turbulent component of the magnetic field Bt. For this discussion 
we therefore define /cad = (^ad ^^^^ that the (isotropic) turbulent power spectrum Tit (k) ~ for 
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Table 1. Results from our fit of Equation (f53l) to the dispersion data for OMC-1. 





Fit Result 






Derived Quantities 




52 (0)b 




Arc 


{B^,)l{BlY {Blf"^ 


(mpc) 




(arcmin~2) 




{^^G) 


16.0 ±0.4 


0.0134 ±0.001 


0.059 ± 0.001 


20.7 


0.28 ±0.01 760 



^Corresponds to the fit result of 7'.'3 ± 0'.'2 at the distance of 450 pc assumed 
for OMC-1. 

^Corresponds to the linear intercept of the broken curves at £ = in the top 
and middle graphs of Figure [2] (see Equation [55]). 

'^Calculated using Equation ([52]l with the SHARP beam radius W = 4'.'7 (or 
FWHM = 11") and A' = 3f5. 

■^Calculated by multiplying the fit result for (0) by N (see Equations [52] 
and [55]). 

'^Calculated using the Chandrasekhar- Fermi equation (see Equation [57]) as- 
suming a density of 10^ cm~^, a mean molecular weight of 2.3, and a velocity 
dispersion of 1.85 km s~^. This estimate is probably not precise to better than 
a factor of a few due to the uncertainty in the density and A' . 
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k > A;ad- It follows that if 6 ^ is some measure of the width of TZt {k) (e.g., its standard deviation) 
then we expect that 5 > (5ad- 



Recent theoretical (ILazarian et al.ll2Q04l ) and observational studies have yielded estimates of 
the order of 1 mpc for 5ad in molecular clouds. Although an observationa l deter mination for the 
ambipolar diffusion scale has yet to be obtained for OMC-1, Li fc Houde (|2008l ) have measured 
5ad ~ 2 mpc for the molecular cloud M17. If we assume for the moment that this value also 
applies to OMC-1, then we find that the turbulent correlation length scale of approximately 16 mpc 
we measured, which is a factor of many larger than the quoted ambipolar diffusion cutoff scale, is 
consistent with the aforementioned expectation that 6 > 6a_d- 



One should always keep in mind, however, that the value for (5ad determined by iLi fc Houde 



(|20Q8l ) pertains to a "Kolmogorov-like" turbulent power spectrum whereas our value for 6 does not. 
That is, given a three-dimensional isotropic turbulent autocorrelation functiorfl TZt (^) from which 
one evaluates the correlation length S, the corresponding power spectrum TZt (k) (of width 6~^) is 
specified by the Fourier transform of TZt (^)- The related Koln iogoroy-like power spectrum is not 
TZt (k), however, but is usually defined as TZk (k) = Airk'^TZt (k) (jFrischlll995l ). It follows that if <5ad 
is the ambipolar diffusion scale pertaining to the Kolmogorov-like spectrum TZk {k) , then we should 
define a turbulent correlation length Sk that also pertains to TZk (k) for a meaningful comparison. 

For example, if we consider the idealization of the Gaussian turbulent autocorrelation function 
of width 6 used for our analysis 



TZt (i) 



(5^ 



then 



TZt{k) = {27Tf^6\~-^'"''" 
and the Kolmogorov-like power spectrum is given by 



(59) 



5/2 x3i.2-y^k^ 



TZk {k) = 2{2TTY''5-'k'e 
Calculating the (square of the) width 5^ of TZk (k) (i.e., its variance) we have 



(60) 



k^TZK (k) dk 
/o°° TZk {k) dk 
= 35-2. 



(61) 



^For this discussion we do not restrict I to the plane of the sky, but allow it to span the three dimensions. 



- 19 - 



Applying this relation to the result obtained for OMC-1, we now find that 5k ~ 9 mpc. This 
value is now a factor of almost two closer to theoretically and observationally expected values for 
(5ad- Furthermore, we should also keep in mind that our assumption of Gaussian functions is 
probably incorrect. The true turbulent autocorrelation function could yield different values for 6 
and (5k, which may be closer to the expected value for (5ad- 

Another fundamental parameter for the characterization of turbulence is the turbulent to large- 
scale magnetic field strength ratio. With the expectation that the magnetic field will be tied to the 
gas through flux-freezing, which should apply for a significant part of the spectrum where k < 5^^, 
this parameter is a measure of the relative amount of turbulent energy contained in the gas. Our 
aforementioned determined value of 0.28 for the square of the turbulent to large-scale magnetic field 
strength ratio seems to indicate that turbulence does not dominate the dynamics in OMC-1. 

As is shown in Equation (f57l) . the square root of this ratio is used to calculate the strength 
plane-of-the-sky component of the large-scale magnetic field with 



1/2 



0.53. 



(62) 



Although our estimate for the plane-of-the-sky component of the large-scale magnetic field of 



1/2 



760 /iG calculated using this ratio cannot be precise to better than a factor of a few 



because of uncertainties in the gas densit y and A^ this value i s reasonable and in line with other 
independent measurements. For example, Crutcher et al. (|l999l ) measured a line-of-sight magnetic 
field strength of 360 /LiG in OMG-1 using CN Zeeman measurements, which probed densities that are 
comparable to those corresponding to our observations. Incidentally, using these values for the two 
components of the magnetic field we can get an estimate of approximately 65° for the inclination 
angle of the large s cale magnetic field to the line of sight. Again, this value is consistent with the 
results obtained by lHoude et al.l (|2004l ) for this object using an independent technique (i.e., 49° in 
the Orion bar and 65° at a location a few arcminutes northeast of Orion KL). 



4.2. Weaknesses of the Technique 

The use of the dispersion function in the polarization angle, while adequately taking into 
account the process of signal integration implicit to dust polarization measurements, allowed us to 
determine some of the fundamental parameters characterizing magnetized turbulence in molecular 
clouds. Our analysis rests, however, on a few assumptions that require some discussion: 

• Our definition for the integrated magnetic field B (r) (see Eq. [10]) does not perfectly mimic 
the measurement process through which dust polarization data are obtained. A consequence 
of this was noted in the discussion that followed Equation ([311) , where it was found that our 
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analysis seemingly introduces a contribution to the angular dispersion that we do not expect 
to be present in actual data. Although a hypothetical analysis that would precisely duplicate 
the measurement process (see Eqs. [39]-[41]) would be desirable, the quest for an analytical 
solution renders this kind of idealization necessary. The same comment could be made for 
the assumptions of isotropy, homogeneity, and stationarity used throughout our calculations. 
Nonetheless, we expect that the results stemming from our treatment of the angular dispersion 
function are successful in the characterization of magnetized turbulence in molecular clouds. 

• We assumed a Gaussian form for the autocorrelation functions characterizing turbulence. This 
assumption is certainly incorrect and as a result it was not possible to obtain a reasonable 
fit when using part of the data where I < 0f2; the shortcomings of our Gaussian model is 
most clearly apparent in that region (see the discussion in Section [3?3]) . We must therefore 
keep in mind that our estimates for the turbulent correlation length 5 and the turbulent to 
large-scale magnetic field strength ratio are correspondingly uncertain to some extent. On the 
other hand our fit is relatively robust, as small changes in the domain used for the fit (i.e., 
the range of values for I) do not lead to significantly different solutions. The same is true if a 
term proportional to £^ is added to the large-scale function to be fitted (i.e., to the right-hand 
side of Eq. [53]). 

• Another source of uncertainty is our modeling of the telescope beam using a Gaussian profile. 
Although the beam size we quote is based on actual measurements taken from chopped mea- 
surements on pointing sources (e.g., Uranus) and is consistent with other SHARP observations, 
we do not possess a detailed map of the telescope beam profile for our set of observations. The 
aforementioned uncertainty is a consequence of the fact that our determination of 5 stems, 
in part, from a comparison of the correlated turbulent component of the dispersion function 
with the assumed Gaussian telescope beam (see the bottom graph of Fig. [2]). The significant 
difference observed between the two, however, makes it unlikely that the ensuing effect on our 
estimate of the turbulent length scale is important. 

• Our analysis has shown that the presence of turbulence in the polarized emission can bring 
about additional angular dispersion to that due to turbulence in the magnetic field. As was 
discussed in Section 13.11 however, our inability to disentangle these two contributions to the 
dispersion has forced us to introduce approximations that essentially brought about the neglect 
of the turbulent polarized emission. It should be noted that this implies that we overestimated 
the amount of turbulence in the magnetic field (measured through (-B^) / {Bq), as obtained 
through our fit to the data), which in turn translates into an underestimate of the large-scale 
magnetic strength when using the Chandrasekhar-Fermi equation. Although this effect is 
probably small, we once again stress that the effect of turbulence in the polarized emission 
could perhaps be advantageously investigated and quantified through simulations. 

It is interesting to note that although we limited ourselves to the determination of only two pa- 
rameters characterizing the turbulent power spectrum, our technique can in principle be used to 



achieve significantly more. We can verify this statement by considering the Fourier transform 6^ (kv) 
associated with the turbulent component (see Eq. [55] . for example), which can be determined from 
Equation (|A12I) 




) 



dk, 



'U 



(63) 



where TZt(kv, ku) is the Fourier transform of Tlt{v,u) = RF,t{v,u) RB,t{v,u). Since the spatial 
frequency component ku of TZt (k^,, k^) is eliminated through the corresponding integration, we can 
equally write 



where TZt (kv) is now interpreted as the two-dimensional turbulent power spectrum. Under the 
assumption of isotropy it would be expected that TZt O^v) is similar in form to TZt {k), where TZt (k) 
is the three-dimensional turbulent power spectrum defined in Section |4?T1 Since the left-hand side of 
Equation ([Mil can be evaluated by taking a Fourier transform of the data (as shown with the symbols 
in the bottom graph of Fig. ^ and that the beam profile H (k^,) is presumably well characterized, it 
follows that the turbulent power spectrum TZt O^v) can readily be determined through the inversion 
of Equation jMJ (using a simple Wiener optimal filter, for example). 

Unfortunately our polarization map does not have enough spatial resolution to allow us to 
perform this analysis (i.e., our beam profile i7(k^) is much too narrow in frequency space), but 
there is no obvious reason why this should not be feasible with higher resolution observations. A 
much better characterization of the turbulent power spectrum would then result. We will address 
this question in a subsequent publication. 



In this paper we expanded our study on the dispersion of polarization angles in molecular 
clouds and showed how the effect of signal integration through the thickness of the cloud as well 
as across the area subtended by the telescope beam inherent to dust continuum measurements can 
be incorporated in our analysis. We correctly accounted for its effect on the measured angular 
dispersion and inferred turbulent to large-scale magnetic field strength ratio. We also showed how 
to evaluate the turbulent magnetic field correlation scale from polarization data and applied our 
results to the molecular cloud OMC-1. For this object, we find a turbulent correlation length of 
5 ~ 16 mpc, a turbulent to large-scale magnetic field strength ratio of approximately 0.5, and a 
plane-of-the-sky large-scale magnetic field strength (-Bq) ^^"^ ~ 760 //G. 

In future papers we will extend our technique to study the possibility for obtaining a detailed 



6' (k.) 



1 



H{k,)fTZtik,), 



(64) 




5. Summary 
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characterization of the turbulent power spectrum, as well as the discussing evidence for anisotropy 
in magnetized turbulence from dust polarization data. 
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A. Derivations 

A.l. The Autocorrelation Function of the Integrated Magnetic Field 

To obtain Equation (fTTI) we start with 



(B (r) • B (r + £)) = / / H {r - a) H {r + £ ~ a') 



-r-TT / / Rf {v,u) Rb{v,u) dz' dz 
^ Jq Jq 



j2 / j2 

a a a a 

(Al) 



from Equations ([TO]) . (fT2]) . (fT3]) . u = \z' — z\, and = |a' — a|. The integral 



1 

I{v)= I I Rp{v,u)R^{v,u)dz'dz (A2) 
^ Jo Jo 

can be transformed into a one-dimensional integral on account of the assumed stationarity of the 
functions present in the integrand (i.e., they are solely a function of \z' — z\, besides v). To accom- 
plish this, after making the change of variable u = \z' — z\, we integrate over the square surface 
delimited by < z < A and < z' < A along the linear path where — A < n < A. This path is 
also perpendicular to a family of strips (along which u is constant) of infinitesimal area 

= A (^1 - du (A3) 

to first order in du. We then find that 

Hv) = ^J^ (l-^)RFiv,u)RB{v,u)du. (A4) 

We can now further make the change of variables a ^ r — a and a' — > r — a', which is then 
inserted with Equation (|A2p into Equation (|Alj) to yield 
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B-B{e)) = jj jj H {a)H {a' + £) (l - Rp iv,u) Rb iv,u) du 



(fa'<fa, (A5) 



where any dependency on r is done away with since u and v are unaffected by this change of variables 
and we integrate over all of space. 

Equation ((24l) can be derived by first noting that 



u 
A 



27r I — I e 



(A6) 



when (5 <C A. This result will arise for all terms containing a turbulent component (i.e., (^t^) 
and/or (F^)) in Equation (lA5]l . 



Setting u = a' + a and v = a' — a we have 



^(iaivK+.r)+iK-ar 



<52 + 2VF2 



+ 



2(,52 + 2VF2)' 



(A7) 



It follows from this relation, the Jacobian related to the coordinate transformation above, and 
Equation ([23]) that 



H (a) H (a' + I) e 



-f^/2(i5^+2W^) 

4(27rW^2)2 
(52 + 2VK2 



(A8) 



which is dependent on only I and not its orientation. Equation (f24l) follows from Equations (fTTI) . 
dll]), (I22|), (1251)- (1281), dMl), and (lA8l). 



A. 2. The Length Scale of the Large-scale Component 



We have used the region where the distance i is less than a few times the beam radius to fit 
the dispersion data to Equation ([43ll . It is therefore appropriate to inquire as to the validity of this 



procedure. To do so, we first note that the autocorrelation function given by Equation (fTT]l can be 
expanded with its Fourier transform as follows (with TZ {v, u) = Rp {v, u) Rb {v, u)) 



(B(r) •B(r + £)) 




(2vr)^ 



(2vr)^ 



H (a) H (a' + £) 

H (a) H (a' + £) 
1 M 



j2 Ij2 

a a a a 



A 



-A 



1 



\u\ 
~A 



(2vr)^ 



du )■ (fa (fa 




H (a) F (a' + i) e^^"<''"'''U^a'd'^ 



1 - ^ ] e^^""dn 



1 - ^ ] e^'^-^dn 



(A9) 



where the Fourier transform of a function is represented by simply replacing the spatial arguments 
by their k-space counterparts (e.g., TZ (v, u) ^ TZ (k^,, ku)), and x = v + uez with v = a' — a and 
u = \z' — z\. But since 



with 



1 

A 



A 

-A 



\u\ 

'a 



AkuU 



du 



1 

A2 
sine" 



rect ( — ^ rect 



A:„A 



U — T 



A 



dT 



'du 



(AlO) 



we can write 



rect (_^) 



1, 
0, 



\t\ < 
\t\ > 



2 ' 
A 
2 ' 



(All) 



(B(r).B(r + £)) 



\H (k,0 



^ / 7^(k^,, /c„) sinc^ 



dku 



e'^^-^d^k^, (A12) 



where \\H = H [K) H* (k^) and sine (x) = sin {x) jx. 
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The effect of signal integration is made explicit in Equation (|A12j) . More precisely, it is seen that 
the integration through the thickness of the cloud heavily filters the spectral components along the 
line of sight to a small set symmetrically located about A;^ = through the presence of sinc^ (ku^) 
and the integration along ku, while the integration across the beam profile on the plane of the sky 
further filters the signal with (k^)|p. For a Gaussian profile as specified by Equation ([23]l we 
have 



^(k.) 



(A13) 



and it becomes clear the only spectral components that remain in the integrated polarization map 
satisfy ky < W~^; the larger the telescope beam the more heavily the signal is filtered spatially. 



It should also be noted that the derivation of Equation (|211) . which was demonstrated in Sec- 
tion lA.ll can also be achieved using Equation (|A12p under the assumption of Gaussian turbulent 
autocorrelation functions (see Equations [21] and [22]) with (5 <C A. 



Restricting Equation (|A12p to its large-scale component, assuming isotropy, and referring to 
Equation ([281) we can write 



1 



2tT roo 



{2Try Jo Jo 



1 



\HikyW 



—j^ \\H{ky)rjoikj) 



1 

2^ 



no ikv,ku)^.^^2 



I 



1 /■ ^0^^^ 2 ( ku^ I 



(A14) 



where TZq {v,u) = Rf,o iv,u) Rb,o {v-,u) and Jq (x) is the Bessel function of the first kind of order 
0. Performing a Taylor expansion about £ = using the identity 



■/o(x)=x;(-i) 



74)" 



n=0 



(n! 



(A15) 



it is clear that, as expected, a Taylor expansion of (a (£)) only contains terms of even power. 
Equations (IA13P and (|A14h show that, once again, the beam filtering process will remove any 
spectral component with ky > in the large-scale dispersion function; its characteristic length 
scale is therefore constrained by, and limited to a few times, the telescope beam radius. Because 
of this beam filtering effect we expect the part of the spectrum at ky < to be dominant in 

Equation (IA14I) . therefore limiting ourselves to a domain where i is less than a few times the beam 
radius (i.e., kyl < 1) to fit the dispersion data with Equation (l43]l implies that 



Jo {kyi) 



(A16) 
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which justifies the method used to model that data and extract the needed information from it. 

B. Data Analysis 

The data for QMC-1 fro m the SHARP polarimeter studied here have been previously published 



by IVaillancourt et alJ (|2008l ). For our purposes we only include data that satisfy the p > 3ap 



criterion, where p is the polarization fraction and ap its uncertainty. 

The angle differences between each and every pair of data points are calculated as 

A^ij = ^i-^j, (Bl) 
and the corresponding distance between each point 



(B2) 



Note that iij = iji so that a map with data points contains only N{N — l)/2 distinct differences. 
Also note that |A$jj| is constrained to be in the range [0,90] degrees. 

These data are divided into separate distance bins with sizes corresponding to integer multiples 
of the grid spacing Ai that results after processing the SHARP map {Ai = 2'.'37; note that a 
SHARP pixel is approximately 4'.'6 and the beam FWHM for our polarization map is approximately 
11"); the bin for 4 (which corresponds to k pixels) covers (4 - A£/2) < £ij < (4 + A£/2). Within 
each bin k we thus calculate the dispersion function with 

1 - (cos iA<^ij))k , for all (4 - Ai/2) < < (4 + Ai/2) . (B3) 

The dispersion function is corrected for measurement uncertainty within each bin according to 
the uncertainty on each A^ij and propagating the measurement uncertainties on both $j and 
available in the data set, as is explained below. However, since it is often the case that iij < W 
some values for and their corresponding uncertainty will be correlated. We have therefore 

used the following relation for the measurement uncertainty 

a\A^,j) ^ a\^i) + a\^j) - 2a{^i)a{^j)e-''^/^'^\ (B4) 

which exhibits the right behavior when iij = and iij > 2W. The data are then corrected for 
measurement uncertainty with 
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(COS (A^>ij))fc n 



(cos(A^^j))fe 
(cos [a(A$,,)])fc 
(cos(A^»i))fc 
l-l(a2(AcI>,,)),' 



(B5) 



where (cos (A<I>jj))^ ^ is the quantity used in our analysis and for the plots in Figure [2l and an even 
probability distribution was assumed for the measurement uncertainty in A$jj. 

Finally, the measurement uncertainties for the dispersion function 1 — (cos (A<I>jj))^ q is deter- 
mined with 



a 



(cos (A$,,)),J = (sin {'^\^^i3))k + 7 (^OS (A$,,))fc {^\^^^3))k ' (^6) 



for all {^k — A£/2) < £jj < (^^ + A^/2). These uncertainties are those plotted in Figure [2], most of 
them are too small to be seen in the figure, especially at the smallest displacements. 



C. List of Symbols 

In this appendix we list and define symbols for important variables and functions appearing in 
the text, as well as give in parentheses the number of the equation where they are first introduced 
or defined. 

• b: normalized magnetic field (Eq. [3]) 

• b: integrated normalized magnetic field (Eq. [3]) 

• correlated turbulent component (Eq. |55| ) 

• B: total magnetic field (Eq. [2]) 

• Bq: large-scale magnetic field (Eq. [2]) 

• Bt: turbulent magnetic field (Eq. [2]) 

• B: integrated total magnetic field (Eq. [TO] ) 

• Bq: integrated large-scale magnetic field (Eq. |10| ) 

• Bt: integrated turbulent magnetic field (Eq. |10| ) 

• e^: unit vector along r on the plane-of-the-sky (Eq. [5]) 
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• e^: unit vector along the line-of-sight within a molecular cloud (Eq. [5]) 

• F: total polarized emission (Eq. [4]) 

• Fq: large-scale polarized emission (Eq. |15| ) 

• Ft: turbulent polarized emission (Eq. |15|)in 

• F: integrated total polarized emission (Eq. [38] ) 

• H: telescope beam profile (Eq. [4]) 

• Jq: Bessel function of the first kind of order (Eq. [A15| ) 

• i: distance between two measurement positions (Eq. [1]) 

• A^: number of independent turbulent cells in the gas column probed by the beam (Eqs. [29] 
and [52]) 

• P: polarized flux (Eq. [41]) 

• P: integrated polarized flux (Eq. [41| ) 

• Q: Stokes parameter (Eq. |39| ) 

• Q: integrated Stokes parameter (Eq. |39| ) 

• r: position vector on the plane-of-the-sky (Eq. [5]) 

• i?B: autocorrelation function of the total magnetic field (Eq. [T2] ) 

• -Rb,o: autocorrelation function of the large-scale magnetic field (Eq. |18| ) 

• i?B,t: autocorrelation function of the turbulent magnetic field (Eq. |18| ) 

• Rp: autocorrelation function of the total polarized emission (Eq. [T3] ) 

• Rf,o'- autocorrelation function of the large-scale polarized emission (Eq. |19| ) 

• Rp^t- autocorrelation function of the turbulent polarized emission (Eq. |19| ) 

• TZ: product of Rp and Rb (in normal space; Eq. [A9| ) 

• TZq: product of Rpfi and Rb,o (in normal space; Eq. [46] ) 

• TZK{k): Kolmogorov-like power spectrum (Eq. |60| ) 

• TZt'. turbulent correlation function (Eq. [58]) 

• U : Stokes parameter (Eq. |40| ) 
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• U : integrated Stokes parameter (Eq. |40| ) 

• W: telescope beam radius (Eq. [23] ) 

• x: three-dimensional position vector (Eq. [5]) 

• z: position along the line-of-sight within a molecular cloud (Eq. [5]) 

• q(£): normalized large-scale function (Eq. |28| ) 

• 5: turbulent correlation length scale of the magnetic field (Eq. |21| ) 

• 5': turbulent correlation length scale of the polarized emission (Eq. [22] ) 

• 5^: Kolomogorov spectral width (Eq. |61| ) 

• A: maximum depth of a molecular cloud along the line-of-sight (Eq. [4]) 

• A': effective depth of the molecular cloud along the line-of-sight (Eq. |34| ) 

• Afcu: spectral width of 7^o (Eq. [47]) 

• A<i> (£): difference in polarization angles between two positions separated by i (Eq. [1]) 

• (•••): average of some quantity (Eq, [l]) 
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Fig. 1. — Normalized autocorrelation function of the i ntegrated polarized flux as calculated usin g 
the previously published 350-^m Hertz data for OMC-1 ([Houde et al.ll2004l : lHildebrand et al.ll2009l ). 
We chose the width at half magnitude to determine the value for the effective depth of OMC-1; we 
therefore have A' ~ 3^5. 
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Fig. 2. — The dispersion function 1 — (cos [A$ [1]]) for OMC-1 using the 350-//m data obtained with 
SHARP. Top: fit of Equation ([431) (sohd curve) to the data (symbols) when plotted as a function of 
the broken curve does not contain the correlated part of the function (see text); middle: same 
as top but plotted as a function of bottom: the turbulent component of the dispersion function 
(symbols), as obtained by subtracting the data points to the broken curve in the middle graph, while 
the broken and solid curves are, respectively, the contribution of the (assumed Gaussian) telescope 
beam alone (i.e, when 6 = 0) and the fit to the data (i.e., with 6 = 7'.'3). 



